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Perturbation series in QCD are generally asymptotic and suffer from so-called infrared 
renormalon ambiguities. In the context of the standard operator product expansion in MS 
these ambiguities are compensated by matrix elements of higher dimension operators, but the 
procedure can be difficult to control due to large numerical cancellations. Explicit subtractions 
for matrix elements and coefficients, depending on a subtraction scale R, can avoid this problem. 
The appropriate choice for R in the Wilson coefficients can widely vary for different processes. 
In this talk we discuss renormalization group evolution with the scale R, and show that it sums 
large logarithms in the difference of processes with widely different R's. We also show that the 
solution of the ^-evolution equations can be used to recover the all order asymptotic form of 
the singularities in the Borel transform of the perturbative series. For the normalization of these 
singularities we obtain a quickly converging sum rule that only needs the known perturbative 
coefficients as an input. This sum rule can be used as a novel test for renormalon ambiguities 
without replying on the large-j3rj approximation. 
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1. The ^-Evolution 

A well known issue in QCD is that perturbative series in a s do not converge, being asymp- 
totic. In perturbative calculations, this feature is caused by contributions from small loop momenta 
that are factorially enhanced due to running coupling effects. Thus an important task in high order 
QCD perturbation theory is to establish to which order the perturbative series can be considered 
"meaningful", or whether it is possible to define subtraction schemes that improve the conver- 
gence properties. A common instrument to quantify the high order behavior of perturbative series 
is the Borel formalism: poles in the Borel transform of the perturbative series in the strong cou- 
pling a s correspond to ambiguities that scale like Aq C d to some integer power p. Predictions 
for physical quantities in perturbative QCD are rendered ambiguity-free in the framework of the 
operator product expansion (OPE) [|I|], sometimes formulated within an effective theory (EFT) con- 
text. Here renormalon ambiguities of Wilson coefficients are compensated by matrix elements of 
higher dimensional operators [Qj. The most common method to determine poles and the corre- 
sponding residues in the Borel transform is based on the computation of massless quark bubble 
graphs inserted in gluon lines and the "naive non-abelianization approximation" [||], also called the 
"large-ray " approximation, which seems to work reasonably well in many examples. 

In an OPE, short and long wavelength contributions are separated into Wilson coefficients that 
are computed perturbatively, and matrix elements of operators that require non-perturbative meth- 
ods. The most common scheme used in the literature for this separation of contributions is the 
MS -scheme. While the MS Wilson coefficients are free from explicit IR divergences to all orders 
in a s , they may still be sensitive to IR physics since they incorporate arbitrarily small loop mo- 
menta, leading to renormalon ambiguities. The resulting numerical cancellations between Wilson 
coefficients and higher dimensional matrix elements are sometimes difficult to control reliably for 
predictions that require small theoretical uncertainties. This problem can be avoided if schemes 
for Wilson coefficients and matrix elements with explicit infrared subtractions are employed. In 
general such subtractions introduce an additional momentum (cutoff) scale, which we denote as R. 
The choice for the scale R depends on the typical distance scale relevant for the quantity described 
in the OPE expansion. It can therefore differ widely for different applications. 

For heavy quark masses such /^-dependent subtractions schemes have been employed for more 
than a decade to avoid the <^(Aqcd) renormalon ambiguity contained in the pole mass [^J and to 
define short-distance mass schemes that can be determined with uncertainties smaller than Aq CD . 
Depending on the process the proper choices for R range from about 1-2 GeV in B -physics or for 
the description of top jets in the resonance region [|], |5]|, to 15 — 20 GeV for top pair production 
at threshold []|], to scales S> m in cases there the MS mass is the proper mass scheme. For OPE 
predictions such subtractions have (apart for the just-mentioned heavy quark masses) only become 
fashionable recently, relevant e.g. for a ambiguity-free definition of the kinetic energy operator X\ 
in the context of HQET [JsJ] , or for an ambiguity-free definition of the gap-parameter that governs 
the first moment of the soft function describing soft radiation in jet production [Q]. 

In this talk we discuss the /^-evolution of matrix elements in schemes with /^-dependent sub- 
tractions. Much like }i, changes in R are related to changes in renormalization scheme. Due to 
a power-like dependence on R, the solutions are substantially different from the common well- 
known logarithmic renormalization group running. We show that the corresponding /?-RGE sums 



2 



The R-evolution ofQCD matrix elements 



Ignazio Scimemi 



large logarithms contained in the difference of two subtractions involving subtractions scales Rq 
and Ri with Rq <C Ri [§]. This resolves e.g. the problem of large logs in the relation of low-scale 
threshold heavy quark masses with R<<^m q and the MS -mass with R k m q . Moreover, taking the 
formal limit R — > in the R-RGE solution, we can recover the all order asymptotic behavior of the 
singularity in the Borel transform that is associated to the subtraction. For the residue of the singu- 
larity we obtain a rapidly convergent sum rule that can be determined from the series coefficients 
of the subtractions without relying on the fermion bubble approximation. This sum rule represents 
a new method to detect renormalons, and even works for renormalons dominated by non-abelian 
effects which are hard to detect with fermion bubbles and naive non-abelianization. In this talk we 
concentrate mostly on the case p=l, e.g. relevant for defining heavy quark masses that are free of 
the ^(Aqcd) renormalon. But we emphasize that our method can be applied in a straightforward 



way for cases with p > 1 as well. For further details we refer to Refs. [|9p and []10fl. 

Consider a low energy non-perturbative matrix element do with mass dimension p, which 
occurs in an OPE for some physical quantity in the MS scheme and has a leading renormalon am- 
biguity of £?(Aq CD ). We can write do as a sum of a matrix element d(R) that is free of the ^(Aq CD ) 
renormalon and a perturbation series (the scheme dependent subtraction) R P 88(R), which contains 
the €?(Aq CD ) renormalon and introduces the subtraction scale R. We have 

/ a (R) \ n 

6o = e(R)+RP86(R), 8d(R) = ^a n (^->) , (1.1) 

where the a n are constants. Since do is /^-independent it is straightforward to determine the R- 
evolution equation for the 0(R), and corresponding /^-anomalous dimension y[a s ,(R)], 

R ik m = ~ RP rt a ^ ' = tj" (^)" ■ (L2) 

Here the coefficients are Jn-i = pci n —1YX=\kak$n-k-\- Since Rd / cIRAq CD = 0, the series for 
7[a v (/?)] is free of the i^(Aq CD ) renormalon. We can extend the situation to the more general 
case where the matrix element is renormalization scale dependent, do = $o(aO- In this case the 
subtraction is also generally renormalization scale dependent, leading to a matrix element 6 (jti , R) 
in an /^-subtraction scheme. Here 6(jj.,R) satisfies both a renormalization group equation in fi, and 



also an /^-evolution equation for \i = R that has as a form identical to eqn.(|1.2|). We note that in 
the /I -dependent situation, it can be quite subtle to construct a subtraction that leads to consistent 
evolution equations in /j. and R. Examples where this has been achieved are the top quark jet mass 
that is relevant for describing the single top invariant mass distribution in top production [§], and 
the gap parameter in the soft function governing dijet production in e + e~ annihilation [J7[J. 

2. General Solution of the R-RGE 



In order to solve eqn. ( |1.2[ ) we need an all order solution of the a r RGE. Consider the a s -RGE, 
R dcC jp = —(xj(R)/{2n) Lr=o A» {a s {R) / An)" , where /3„ are /3-function coefficients in an arbitrary 
scheme for the strong coupling. A solution of this equation is given by 

ln ^ = f = f°dtb(t) = G(to)-G( h ), (2.1) 

Rq J ao P[CCr\ Jh 
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where a,- = a 4 . (/?,•), a R = a s (R), t t = —2n/(po(Xi), and t = -27r/(fta s ). Here 

G(t)=t+b l ln{- 



'2 »3 
+ — + ... 



b3_ 
2t 2 



(2.2) 



3 = (j3 1 3 -2A)J3iJ3 2 + 
t) . From eqn. (Oh 



t t 2 t 3 ' ^ ' t 

where for the first few orders S a = ft /(2ft 2 ), £ 2 = (j3 2 -ftft)/(4ft 4 ), and b 
J3 ( 2 J83) / (8j8q ). For later convenience we also define G 2 (f) = G(/) 
one immediately notes that, 7? exp G(t) = Ro exp G(fo) = Aqcd- This equality demands Re G ^ to be 
a constant of mass dimension one, which we have defined to be Aqcd and one can easily check that 
this definition corresponds to the familiar definition of Aq CD at N*LL order. The solution above is 
valid for an arbitrary mass-independent scheme for the j3 -function, although usually one uses the 
MS scheme where the j3, are known to four-loops [12]. Here we present the solution for the /?-RGE 
of eqn. (1.2) for p = 1. For arbitrary p, see [PH. Integrating equation (L2) gives, 



e(Ri)-e(Ro) 



Ro 



dRy[a s (R) 



dtb(t) 



Aq C d^ GW 



y[t] 



(2.3) 



where we have used the solution of the a $ -RGE and defined y[t] = y[a s (R)] using the substitution 
t = —2k/ (fta s (/?)). In order to write the solution in a closed form we define 



b(t)e- G2{t) y[t] 

72 



7=0 



s 2 



(S 1 +S 2 )-^+[(i+S 1 )S 2 + (S|+S 3 )/2] 



, where So 



7b 
2ft' 



(2ft 



(2ft) 
7o 



(bi+b 2 ) 



2ft 



(2ft) 



Now substituting eqn. ( |2.4[ ) into eqn. (| 

tion with the standard convention of integration above the cut, f t °°dt' ( 
we get the N^LL solution that is given by 



(2ft) 

|) and using the definition of the incomplete gamma func- 



t '\c-l„-t' 



T(c,0, 



[e{Ry)-e{R Q )\ 



N LL 



(A^^Sji-iye^m-j-hA) 



n-j-b h t Q )] ■ (2.5) 



Here Aq CD is the N^LL solution of: R expG(f) = Aqcd - The integrand in eqn. (2.3) has an essential 
singularity and a branch point at t = due to infinitely many negative powers of t and a branch 
cut on the positive real axis due to the (— t)~ b[ . For Rq,Ri > Aqcd we have to,t\ < 0, so the range 
of integration is away from the branch cut and the integral is well denned. To write a closed form 
expression we have split the integral into a difference of two incomplete gamma functions, and with 
the e ,nbl factor this solution is real. Note that only the difference of the gamma functions appearing 
in eq. (2.5) is independent of the convention for treating the cut. 



To see what kind of perturbative terms the solution in eq. ( |2.5| ) contains, we take the k = 
case (LL) and recall the asymptotic expansion for t — > —00, T[c,t] a = m e~ l t c ~ l Y^=o (^yr^i-c) • 



For c = this yields A^rfO,?] "= -2fl£~ =0 2"n! f p a s (R) / (4n) J , where we used the LL 



asym 



relation Ag CD e~ 



R. This is a divergent series, but for the LL solution of the RGE 



[e{R l )-e{R Q )\ 



LL 



-TfoRi y r ft«i 

2ft ^ L 2% 

y Ro y rftai 

" 2ft [ 2ti 



n+l 



Ro 



n+l 



I 

k=n+\ 



k=0 

k\ Ro 



Ri 
Ro 



(2.6) 
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which is convergent since f5oCC s (Ri)ln(Ri / Ro) / (2n) < 1. Eqn. (2.6) displays the problem of large 
logs in fixed order perturbation theory for R\ 3> Rq. The RGE in R encodes IR physics from 
the large order behavior of perturbation theory It causes a rearrangement of the IR fluctuations 
included in 6(R) in going from R to Ri without reintroducing a renormalon. 



3. The R-RGE and the Leading Renormalon Ambiguity 



In this section we demonstrate a novel application of the i?-RGE by deriving the leading renor- 



malon ambiguity of do of eqn. (1.1) in the Borel plane and providing a sum rule to perturbatively 
determine the normalization of the renormalon ambiguity. A sum rule of this kind for u = 1/2 
ambiguity of the pole mass was first provided in our earlier work [^]. 

We note 1 that in the complex 7?-plane away from the positive real axis, limR_>o OC s (R) = 0. 
Therefore 8d(R) of eqn. ( pTT| ) vanishes at R = and we get lims^o^(^) = 6b- It i s important 
to note that though a s (R = 0) is single- valued, but the process of taking the limit 7? — > is not 
unambiguous due to the branch cut on < R < Aqcd- It is the limiting procedure that reintroduces 
the ambiguity. Thus, taking limit Rq — > in eqn. ( |Q| ) we get 



e(Ri)-e Q 



-Aqcd £ Sj 
;'=o 



dt- 



(3.1) 



The LHS has a renormalon ambiguity in do, which on the RHS is in the integral, the integrand has 
a branch cut on the positive real axis in the complex ? -plane, and the integration must be performed 
by either going above the positive real axis or below it. Using the asymptotic expansion for this 
integral, which is the same as that for T(c,t), and expanding e Gl ^ = Yd=o tz^i > we g et a power 
series expansion in l/(— /) = poCC s {R)/(2n), 

6(R)-6o = -R Y S j8t m+J+n + h) 



On taking the Borel transform 2 , after some algebra, we get 



B(u) = 2R 



1=0 



r(i+&i-i) 



where 



1/2 



limP 1/2 



N k LL 



i im y 



Sj 



+bi+j) 



(3.2) 



(3.3) 



and Q(u) is a function convergent in a finite neighborhood of u = 1/2. The difference, 6(R) — 
O , can be obtained by the inverse Borel transform, d(R) - % = £du B(u) e~ u 43r /(ft« s (K)) . This 



1 Eqn. ( 1. 1 ) has the solution that as / goes to infinity, R — > with a phase e ±mhl . That is to say R cannot approach 
from the positive real axis. R must go to either from above the positive real axis or below it. This is so because there is 
a cut on the positive real axis for < R < Aqcd- This cut corresponds to the branch cut in the f -plane for the integrand 
of the equation ( |2.3| ). Nevertheless, t as a function of R is single- valued at R = 0, corresponding to a s (R = 0) = 0. 

2 Taking the Borel transform amounts to making the replacement \/t n+l — > 2(2«)"/r(n+ 1). In our convention, 
which is also the most widely used, u is Borel conjugate of fta.,/ (4s). 
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integral is ambiguous because eqn. (3.2) has a branch cut for u > 1/2. Thus we have obtained, 



in eqn. (p.2|), the structure of the leading renormalon ambiguity in 6q, represented by the second 



term on the RHS, along with its normalization Pin- Here P 1 / 2 is an absolutely convergent infinite 



series [|9J] constructed from the coefficients Sj that appear in the solution of the R-RGE (eqn. ( |2.4[ )). 



The Sj are in turn determined by the /^-anomalous dimension coefficients y„ (eqn. (|1.2[)). Therefore 
is perturbatively determinable in terms of the original coefficients a n of 8d(R). 
For any given d(R) the process described above can be repeated with a rescaling R — > XR. 
This does not change the normalization of the renormalon ambiguity, since it amounts to only 
changing the scale R, or equivalently changing the scheme. In this new A -dependent scheme, we 
have A -dependent Sj, where the first few are given by 

iSo(A) = ASo , Si (A) = A>Si — 



2A> 

Alog(A)(-a 2 + ai(%+2)j3Q + aij3olog(A)) 
2$ 



S 2 (A) = XS 2 + -6^A-»2Tu U , 2 -r^Tu 1Ml ^^, (3 4) 



As a result P]y 2 LL depends on A, but Pi/ 2 will be A -independent. This gives a way to assign errors 

and test for convergence of the partial sum ^J/ 2 LL . Distinguishing P 1( / 2 ^ and Pi/ 2 = provides 
a test for the presence of a renormalon ambiguity. We propose to keep 0.5 < A < 2, so that the 
perturbative expansion is well behaved. 

4. Applications of the Renormalon Sum Rule 

As an application of the sum rule, we calculate the residue of the renormalon ambiguity of 



the heavy quark pole mass. Using eqn. ( |3.3[ ) along with Sj's from eqn. (p.4|), we obtain Pj/ 2 as a 
function of scaling A in MSR [[| |Tl|,|lO[], static [g, |l(j] and PS [{3|] mass schemes up to NNLL order. 
This is shown in Figure [j], where the band that envelops P Y i 2 curves for these schemes (with rif = 5 
light flavors) is plotted at LL, NLL and NNLL order. Since the pole mass renormalon ambiguity is 
scheme independent, we expect the sum rule to converge to a scheme independent number. This is 
clearly seen from the figure, since the band at the NNLL order is significantly narrower with weak 
A-dependence. From the narrow red (NNLL) band we estimate that P™"™ = 0.47 ±0.10, which 
is the normalization of the order Aqcd renormalon ambiguity in the heavy quark pole mass. For 
comparison, the bubble chain approximation gives Pfff s = 0.80, overestimating by a factor of 2. 
Our values are consistent with the determination of the normalization of the renormalon ambiguity 



in Ref. [jlj. 

The Wilson coefficient C cm (niQ , pL ) of the chromomagnetic HQET operator h v o ■ Gh v , is known 
to have a ^(Aqcd/wiq) renormalon. We will test for this renormalon using our renormalon sum 



rule and determine its size. C cm is known to three loops fll3[ ] and at ju = rag, it is given by, 
Ccm(niQ,niQ) — 1 = Y°°=i a i ( "'4^ ) ■ W e are looking for n\(2po) n growth (p = 1) in the coef- 



ficients of a s (/?)"(mg). To test this, we take the known MS scheme results ai = 8.66703, a 2 = 
(350.347 - 30.6037?!/), a 3 = (21985.2 - 3470.72n/ + m.Snj) [0] and substitute in eqn. Q, 
which we then use in eqn. ( J3.3| ) to get the P X j 2 sum rule as a function of A. This sum rule is 
plotted in Figure [T] for rif = 3 (thin lines) and rif = 4 (thick lines) light flavors at the LL(dotted), 
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pmass 
r l/2 



NNLL 0i 



0.6 




nchromo 
1.0 




0.4 



0.2 



1 LL 



NNLL, 

—4^ 



NLL 



\ 2.0 0.5 1.0 1.5 2.0 2.5 \ 3.0 

Figure 1: The Pjm sum rule for the heavy quark pole mass (left) and the chromomagnetic operator Wilson 
coefficient (right). 



NLL(dashed) and NNLL(solid) orders; convergence is quite evident and from the A dependence 



0.72 ±0.09 (n/ = 3), and p^ mo = 0.71 ±0.07 (n f = 4). The sum rule 



we estimate p^ omo 

yields a clear renormalon. For further details and discussions regarding concepts presented in this 
proceeding, including some direct applications of the /?-RGE, see [[| [Tc| ]. 
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